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Abstract 

A new approximation scheme to the centrifugal term is proposed to obtain the I ^ bound-state 
solutions of the Schrodinger equation for an exponential-type potential in the framework of the 
hypergeometric method. The corresponding normalized wave functions are also found in terms 
of the Jacobi polynomials. To show the accuracy of the new proposed approximation scheme, 
we calculate the energy eigenvalues numerically for arbitrary quantum numbers n and I with two 
different values of the potential parameter ctq . Our numerical results are of high accuracy like the 
other numerical results obtained by using program based on a numerical integration procedure for 
short-range and long-range potentials. The energy bound-state solutions for the s-wave (/ = 0) 
and (To = 1 cases are given. 
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I. INTRODUCTION 



The exact analytic solutions of the wave equations (nonrelativistic and relativistic) are 
only possible for certain potentials of physical interest under consideration since they contain 
all the necessary information on the quantum system. It is well known that the exact 
solutions of these wave equations are only possible in a few simple cases such as the Coulomb, 
the harmonic oscillator, pseudoharmonic potentials and others [1-5]. Recently, the analytic 
exact solutions of the wave equation with some exponential-type potentials are impossible 
hr I ^ states. Approximation methods have to be used to deal with the centrifugal term 
like the Pekeris approximation [6-8] and the approximated scheme suggested by Greene 
and Aldrich [9]. Some of these exponential-type potentials include the Morse potential 
[10], the Hulthen potential [11], the Poschl- Teller [12], the Woods-Saxon potential [13], the 
Kratzer-type and pseudoharmonic potentials [14], the Rosen-Morse-type potentials [15], the 
Manning- Rosen potential [15-22], other multiparameter exponential-type potentials [23,24] 
and hyperbohcal potential [25-27]. 

In this work, we attempt to study another exponential-type potential called the hyper- 
bohcal potential [25-27] 



where D, a and ctq are three positive parameters. It is indicated in [25] that this exponential- 
type potential is closely related to the Morse, the Kratzer, the Coulomb, the harmonic 
oscillator and other potential functions. The properties and applications of this potential 
are given in [25,26]. It is known that for this potential the Schrodinger equation (SE) can 
be solved for the s-wave, angular momentum quantum number / = 0. However, for a general 
solution, it is need to include some approximations if one wants to obtain analytical or 
semianalytical solutions to the SE. For the I ^ case, the potential (1) can not be solved 
exactly without an approximation to the centrifugal term [22]. Hence, in the previous 
papers, several approximations have been developed to find better analytical formulas for 
the hyperbohcal potential [27]. 

Our aim in this work is to attempt to study the arbitrary /-state solutions of the 
Schrodinger equation for the hyperbolical potential. In order to improve the accuracy of 
our previous approximation [20,21], we propose and apply a new approximation scheme for 




(1) 
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the centrifugal term in the form: 

1 



„ - 4q; 



1 - e-^"'- V 1 - 



(2) 



where Cq is a proper shift to be found by the expansion procedures. Thus, with this new 
approximation scheme, we calculate the I ^ energy levels and wavefunctions of the hyper- 
bolical potential using the hypergeometric approach (Nikiforov and Uvarov) N-U method. 
This method has shown its power in calculating the exact energy levels for some solvable 
quantum systems [13,14,20-22]. The approximation given by (2) has proved its power and 
accuracy over the other currently used approximations in literature [16-21]. It has been 
applied recently on the Manning-Rosen potential [22] and has also proved its power and 
efficiency when compared with the other numerical simulations for the non- approximated 
problem used to calculate the energy bound states. It provides good results which are in 
agreement with the numerical integration method by Lucha and Schoberl [28]. 

The paper is organized as follows: In Section II we breifiy present the N-U method. In 
Section III, we present the new proposed approximation scheme and apply it to calculate 
the Z-wave eigensolutions of the SE for the hyperbohcal potential by the N-U method. In 
Section IV, we present our numerical results for energy eigenvalues numerically for arbitrary 
quantum numbers n and / with two different values of the potential parameter (Tq. Section V, 
is devoted to study two special cases, the s-wavc (/ = 0) case and the ctq = 1 exponential-type 
potential. Finally, we make a few concluding remarks in Section VI. 



II. THE NIKIFORV-UVAROV METHOD 

The Nikiforov-Uvarov (N-U) method is based on solving the second-order hnear differen- 
tial equation by reducing it to a generalized equation of hypergeometric type [29] . In this 
method after employing an appropriate coordinate transformation z = z{r), the Schrodinger 
equation can be written in the following form: 

where (t{z) and '^{z) are the polynomials with at most of second-degree, and t{s) is a first- 
degree polynomial. The special orthogonal polynomials [29] reduce Eq. (3) to a simple 
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form by employing ipniz) — (f)n{z)yniz) , and choosing an appropriate function (pniz). Con- 
sequently, Eq. (3) can be reduced into an equation of the following hypergeometric type: 

a{z)y';{z) + r{z)y'^{z) + Xy^iz) = 0, (4) 

where t{z) = t{z) + 2tt{z) (its derivative must be negative) and A is a constant given in the 
form 

A = A„ = -nT'(^)-^^i^^a"(^), n = 0,l,2,... (5) 

It is worthwhile to note that A or A^ are obtained from a particular solution of the form 
y{z) = yn{z) which is a polynomial of degree n. Further, yn{z) is the hypergeometric-type 
function whose polynomial solutions are given by Rodrigues relation 

where is the normalization constant and the weight function p{z) must satisfy the dif- 
ferential equation: [29] 



In order to determine the weight function given in Eq. (7), we must obtain the following 
polynomial: 

a'{z) - r{z) ^ lfa'{z)-T{zy ' 



^(^)- ' \ V 2 ) + (8) 

In principle, the expression under the square root sign in Eq. (8) can be arranged as the 
square of a polynomial. This is possible only if its discriminant is zero. In this case, an 
equation for k is obtained. After solving this equation, the obtained values of k are included 
in the N-U method and here there is a relationship between A and A; by A; = A — 7r'(2;). After 
this point an appropriate (f)„{z) can be extracted from the differential equation: 

<P'i^) - (^) H^) - 0. (9) 



a(z) 

III. ANALYTICAL SOLUTIONS 

A. An Impoved Approximation Scheme 



The approximation is based on the expansion of the centrifugal term in a series of expo- 
nentials depending on the intermolecular distance r and keeping terms up to second order. 
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Therefore, instead of using the approximation in [9,11,19], we use this choice of approxima- 
tion: 



1 



-2ar 



— ^ {2af [co + v{r) + v^{r)] , v{r) = - 



-2or ' 



17 ^ (2«)' 



Co + 
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:,2ar 



+ 



32ar 



(10) 



which has a similar form of the hyperbohcal potential. Changing the coordinate to x by 
using X — {r — ro)/ro, one obtains 



{l + x)-'^j' 



Co + 



1 



+ 



1 



, 7 = 2aro 



(11) 



e7(l+^) _ 1 ' (g7(l+x) _ 

and expanding Eq. (11) around r = ro (x = 0), we obtain the following Taylor's expansion: 

l-2x + 0{x^)^j^ (^co + 
1 3 
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1 1 

e^"^ ^ (e-r - 1) 

2 



eT - 1 (e7 - 1)^ (e7 - 1) 



x + 0{x^), 



(12) 



from which we obtain 
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1 1 
Co + ^ r + 



eT - 1 (eT - 1)^ 



= 1, 



7 



+ 



3 2 

+ 



2. 



(13) 



,eT-l ■ (e^-l)^ (e-^- 1)-^ 
Therefore the shifting paramete, cq, can be found from the solution of the above two equa- 
tions as: 

Co = A 7 — ^—^ = 0.0823058167837972, (14) 

72 eT - 1 (eT - 1)2 ^ ^ 

where e is the base of the natural logarithms, e = 2.718281828459045 and the parameter 

7 = 0.4990429999. Hence, we have the following substitution for the centrifugal term: 



lim4Q; 



a-»0 



1 a 
+ 



72 - 1 (eT - 1)2 1 - e 



-2ar 



1 - 6-2°'- 



(15) 



Finally, it is worth to note that, in the case if Cq = 0, the approximation given in Eq. (10) 
is identical to the commonly used approximation in the previous works [9,11,19-21]. 
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B. Energy Eigenvalues and Eigenfunctions Solution 



To study any quantum physical system characterized by the empirical potential given in 
Eq. (1), we solve the original SE which is given in the well known textbooks [1,2] 

(l^ + ^W) ^(r'^' ^) = ^^(r'^' '^)' (16) 

where the potential V{r) is taken as the hyperbolical potential (1). Using the separation 
method with the wavefunction ip{r,9,(f)) — r~^R{r)Yi,m{9,4>), we obtain the following radial 
Schrodinger eqauation as 

1 - (Jo 



2 



- ^a%l + 1) Co + 



\2 



(1 - e-'^'^'^y 

X R^i{t) = 0. (17) 

Since the SE with the above hyperbohcal potential has no analytical solution for Z-waves, 
we have used the approximation to the centrifugal term given by case 1. The other approx- 
imations will be left for future investigations. To solve it by the N-U method, we need to 
recast Eq. (17) into the form of Eq. (3) changing the variables r ^ z through the mapping 
function r — f{z) and making the following definitions: 



z = e-2"^ u = s' = W-f^ + AEi, E^i < AEi, AEi = l{l + l)co, (18) 



we obtain the following hypergeometric equation: 

d'^R{z) (1 - z) dR{z) 1 
dz^ ^ z{l-z) dz ^ [^(1 - z)f 

X {-s'^ aof + [2u (1 - (Jo') + 2s''' - 1(1 + l)]z- [u (1 + aof + s'^] z'} R{z) = 0. 

(19) 

It is shown from Eq. (18) that for bound state (real) solutions, we require: 

0, when r — > oo, , , 

(20) 

1, when r — > 0, 

and as a result the radial wavefunctions Rni{z) for the values of z given in Eq. (18). To 
apply the N-U method, we compare Eq. (19) with Eq. (3) and obtain the following values 
for the parameters: 

t{z) — 1 — z, a{z) — z — z'^, 



a{z) + aof + e'^] + \lv (l - a^) + 2£'2 _ + 1)] ^ - s'^ _ (i _ ^o)^ . (21) 

If one inserts these values of parameters into Eq. (8), with (j'{z^ = 1 — Iz^ the following 
linear function is obtained 

7r(2;) = -- ± -A/a^I^TTo^TTa^, (22) 

where 02 = 1 + 4 [e'^ + (1 + ctq)^ - A;] , ai = 4 {A; + i(i + 1) - 2i/ (1 - cr^) - 2e'^} and Oq = 
4 [s'^ + i/ (1 — (7o)^] . According to this method the expression in the square root has to be 
set equal to zero, that is, A = a2Z^ + a\z + ao = 0. Thus the constant k can be determined 
as 

k=-\ [(1 + 2/)' + 16z/ao(ao - 1) - l] ± + 25), (23) 

where 

/3 = £'2 + ^(i_^^)2^ ^ ^ 1 _^ ^ ^l6i/a2 + (l + 20^ 
In this regard, we can find four possible functions for 7r(2;) as 



(24) 



X 



(2/3-25-1)^-2/3, for A; = -[(l + 202 + 16i/(7o((7o-l)-l]/4 + /3(l + 25), 

(2/3 + 25 + 1)^-2/3; for A; = - [(1 + 2/)^ + 16i/(7o((7o - 1) - 1] /4 - /3(1 + 25). 

(25) 

We must select 

A; = -|-^[(2/3 + 25+l)z-2/3], (26) 
in order to obtain the polynomial, t{z) — t{z) + 2ti{z) having negative derivative as 

t{z) = 1 - 2z - [2/3 + 25 + 1] z - 2/3, t\z) = -(2/3 + 25 + 3). (27) 

We can also write the values of X — k + n'{z) and A„ = —nT'{z) — '^^'^~^'> a"{z), n = 0, 1, 2, ... 
as 

A = - J [(1 + 2lf + IQuaoiao - 1) - l] - [/3 + 5 + 1] , (28) 

A„ = n(n + 2/3 + 25 + 2), n = 0, 1, 2, • • • (29) 

respectively. Additionally, using the definition of A = A„ and solving the resulting equation 
for allows one to obtain 



n + l)2 + /(/ + l) + (2n + l)5-4i/ao(l-ao) / ^iEnl , . ^ , .2 

2(n + 5 + l) ^V"^ + '^^' + "^'""°^- 

(30) 



Hence, we obtain analytically the following discrete bound-energy levels 



Eni = (1 - (JofD + 



2a^n%l + 1) 



72 - 1 (e^ - ly 



2*2 



{n + 1)2 + l{l + 1) + (2n + 1)5 - 4i/(Jo(l - (Jq) 



H 2 



,0<n,Z<oo (31) 

^i I 2(n + (5 + l) J 

where n — 0,1,2, ■ ■ ■ and Z signify the usual radial and angular momentum quantum num- 
bers, respectively. 

Let us now find the corresponding radial part of the normalized wave functions. Using 
a{z) and 7r{z) in Eqs. (21) and (26), we obtain 



p{z)^z''f'{l-zf'+\ 



(32) 
(33) 



(34) 



The functions yni{z), up to a numerical factor, are in the form of Jacobi polynomials, i.e., 
yni{z) ~ Pi^^'^^^^\l - 2z), valid physically in the interval (0<r<oo^O<z<l) [30]. 
Therefore, the radial part of the wave functions can be found by substituting Eqs. (32) and 
(34) into Rni{z) = 4>{z)yni{z) as 



Rni{r) = iV,;e-2"'^'-(l - e-2-)i+^p^2A25+i)(i _ 



(35) 



where /3 and S are given in Eq. (24) and A^^; is a normalization factor to be determined from 
the normalization conditionrwhich gives [20-22] 

1 



s(n) 



X 



E 

p,r=0 



_ (-1)" r{n + 25 + 2)r{n + 2/3 + 1)^ 
~ r{n + 2(3 + 25 + 2) 

l)P+'T{n + 2p + r-p+l){p + 25 + 2) 



p\r\{n - p)!(n - r)\T{n + 2(3 -p + 1)1(2/3 + r + l){n + 2(3 + r + 25 + 2)' 



(36) 
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IV. NUMERICAL RESULTS 



To show the accuracy of the new approximation scheme, we calculate the energy eigen- 
values for various n and I quantum numbers with two different values of the parameter (Jq. 
The results calculated by Eq. (31) are compared with those obtained by a MATHEMAT- 
ICA package programmed by Lucha and Schoberl [28] as shown in Table 1 for short-range 
potential (small a) and long-range potential (large a). It provides that the new proposed 
approximation scheme to the centrifugal term in Eq. (2), even when the potential parame- 
ter a becomes large, produces energy eigenvalues of high accuracy like the other numerical 
methods [28] . Consequently, this is also an illustration to assess the validity and usefulness 
of our present approximation. Further, it is quite simple, computationally efficient, reliable 
and accurate. 



V. DISCUSSIONS 



We have used the hypergeometric method (N-U) to solve the radial SE with the 
exponential-type potentials for arbitrary Z-states. We have derived the binding energy spec- 
tra in Eq. (31) and their corresponding normalized wave functions in Eq. (35). 

Firstly, let us attempt to study the s-wave case (/ = 0). Hence, the energy eigenvalue and 
the radial eigen function solutions; (31) and (35), reduce to the following forms: 



2fc2 



En = D{1 - aof - 



1^ 



{n + ly + (2n + l)Si - 4uao{l - ctq) 
2{n + Si + l) 



, < n < oo, (37) 



and 



Rni{r) = Nnie-^^'^^'il - e-2«'-)i+5ipm,25i+i)(i _ 2e-2«'-), 



respectively, with 



-1 + 11 + 



a 



(38) 



(39) 



and Nni is given in Eq. (36). This result is in agreement with Ref. [27]. 

Secondly, we further discuss another special case ctq = 1. As a result, the potential (1) 
reduces to 



V{r) 



— 4ar 



ADe 
(1 - e-2«'^) 



2 ■ 



(40) 
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The corresponding energy levels and radial wave functions are given by 



2a^h'^l{l + 1) 



nl 



and 



72 eT - 1 (e-y - ly 



2fe2 



1^ 



{n + If + l{l + 1) + {2n + 1)52 



2{n + 52 + 1) 



1) 



Rni{r) = N^ie-^'^^^'il - e-2"'-)i+'52pP/?2-252+i)(i _ 2e-^'''), 



(42) 



respectively, with 

/?2 = 



llEnl , /(/ + 1) 1(1 + 1) 1(1 + 1) ^ 1 

H ^ 7 : — —, 02 - 



2rf ' 72 e^-1 (e^-l)2' 2 



-1 + 



+ (l + 2/)2 



(43) 

and Nni is given in Eq. (36). This result essentially coincides with that of Manning- Rosen 
potential with the special case A = and D — "^(^~")^ shown in our recent work [22]. 



VI. COCLUDING REMARKS 



The arbirary /-wave solutions of the SE with an exponential-type potential have been 
obtained approximately by proposing an improved shifted approximation to the centrifugal 
term. It is found that the normalized wave functions can be expressed by means of the 
Jacobi polynomials. With this approximation scheme, we can easily build an analytic for- 
mulations [Eqs. (31) and (35)] and use it to evaluate eigenvalues and eigenfunctions. This 
approximation scheme has also been used with great success in problems which do not have 
exact solutions for / 7^ case with exponential-type potentials like Manning-Rosen poten- 
tial and hyperbolic potential. Essentially, two special cases for the s-wave case (/ = 0) and 
(Tq — 1 and found that these results have been reduced to those given in [25-27]. To show the 
accuracy of our results, we have calculated the eigenvalues numerically for arbitrary n and 
I with two different values of the parameter uq. We found that the results obtained by (31) 
are in good agreement with those obtained by using the MATHEMATICA program based 
on the numerical integration procedure for short-range potential (small a) and long-range 
potential (large a) [28]. As a demonstration of the accuracy of our results, Table 1 shows 
that the estimated energy eigenvalues can be computed up to 0.001 — 0.13 %. However, 
the accuracy of the energy eigenvalues in [27] is computed up to 0.051 — 1.0 %. Therefore, 
the accuracy of the present model reaches up to 10 — 50 times better than the estimations 
provided by [27]. 
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TABLE I: Energy eigenvalues as a function of the parameter a for 2p, 3p, 3d, 4p, 4(i, 4/, 5^», 5ci, 5/, 
hg^ 6p, 6(i, 6/ and states for (Tq = 0.1, 0.2 and D = 10 in atomic units {h = jj, = 1). 
(70 = 0.1 (70 = 0.2 

states a present Dong et al [27] Lucha et al [28] present Dong et al [27] Lucha et al [28] 
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